Absence of vertex correction for the spin Hall effect in p-type semiconductors 
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We calculate an effect of spinless impurities on the spin Hall effect of the Luttinger model rep- 
resenting p-type semiconductors. The self-energy in the Born approximation becomes diagonal in 
the helicity basis and its value is independent of the wavenumber or helicity. The vertex correction 
in the ladder approximation vanishes identically, in sharp contrast with the Rashba model. This 
implies that in the clean limit the spin Hall conductivity reproduces the value of the intrinsic spin 
Hall conductivity calculated in earlier papers. 
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Understanding the dynamics of spins is a long-standing 
subject in semiconductor physics. Though there have 
been lots of experimental and theoretical studies toward 
this goal, there remain many obstacles to overcome before 
we can manipulate spins at our disposal in semiconduc- 
tor spintronics devices^. One of the obstacles is efficient 
spin injection into semiconductors, and there have been 
many atempts for it. For example, spin injection from 
ferromagnetic metals has low efficiency, because of the 
band structure mismatch between semiconductors and 
metals, in particular by the conductivity mismatch 3 ". On 
the other hand, while spin injection from ferromagnetic 
semiconductors such as (Ga,Mn)As^ provides relatively 
high efficiency^, the Curie temperatures of such ferro- 
magnetic semiconductors still remain lower than room 
temperature. Thus the spin injection, which constitutes 
the very first step for spintronic devices, is still under 
intensive investigation. 

In such circumstances, recent theoretical proposals of 
an intrinsic spin Hall effecliSii have attracted much at- 
tention. In these works, it is theoretically predicted 
that an external electric field induces a dissipationless 
spin current in semiconductors. This intrinsic spin Hall 
effect is predicted for the Luttinger model represent- 
ing bulk p-type semiconductors^, and for the Rashba 
model representing two-dimensional n-type semiconduc- 
tors in heterostructuret This intrinsic spin Hall ef- 
fect is caused by the Berry phased in the momentum 
space. Berry-phase structure in momentum space is 
comprised in the electronic Bloch bands&iSiii, and en- 
dows the carriers with an anomalous velocity. Due to 
the spin-orbit coupling, the carriers will then have spin- 
dependent trajectories in the presence of an external elec- 
tric field, and it amounts to the intrinsic spin Hall ef- 
fect. In bulk p-type semiconductors, the valence band 
has a spin splitting, giving rise to spin Hall effect by hole- 
doping as proposed by the author and his collaborators 6 , 
whereas in bulk n-type semiconductors, the conduction 
band is doubly degenerate (if the Dresselhaus term is 
neglected), and the spin Hall effect does not emerge. In- 
stead, by making the heterostructure from the n-type 
semiconductor, the Rashba spin splitting 12 *! 3 " arises, and 
it gives rise to the spin Hall effect proposed by Sinova 
et ali Although an unambiguous experimental detec- 



tion remains to be done, many theoretical studies have 
appeared since the proposals of the intrinsic spin Hall 
effect 14 ! 15 ! 16 !^! 18 ! 13 ! 2 "! 21 ! 22 ! 23 ! 24 ! 25 ! 26 . We note that these 
are different from an extrinsic spin Hall effect caused by 
spin-dependent scattering by impurities, which might be 
assumed as a main source of spin Hall effect in earlier 
paper a 2 Ti 2 ?i 2 ? . 

In order to argue possibilities of experimental detec- 
tion, theoretical investigation of disorder effects is highly 
desired. In this paper, we mainly focus on the Luttinger 
model representing p-type semiconductors^, and study 
disorder effects on the spin Hall effect. We shall also dis- 
cuss on the Rashba model. For the Luttinger model, 
the effect of self-energy broadening was calculated by 
Schliemann and Lossii, and no other papers on disor- 
der effect have appeared at present. In this paper, we 
study the effect of randomly distributed spinless impuri- 
ties with short-ranged potential, by calculating both the 
self-energy correction in the Born approximation and the 
vertex correction in the ladder approximation. The cal- 
culation itself is analogous to that by Inoue et al. for 
the Rashba modelS 4 *^!. The resulting self-energy is of a 
simple form, independent of the wavevector k. Remark- 
ably, the vertex correction turns out to vanish identi- 
cally, even away from the clean limit. Thus the spin Hall 
conductivity becomes identical with that of the intrinsic 
one calculated in ref. ||J We note that it is in sharp con- 
trast with the Rashba model, where the vertex correction 
completely kills the spin Hall effect within the Born and 
ladder approximations in the clean limit 24-31 . 

The Luttinger Hamiltonian 30 is written as 

H ° {k) = £ ( (71 + \ l2)e ~ 272(k ' S)2 ) (1) 

where S = (S x , S y , S z ) are the spin-3/2 matrices, 
k = [k x , k y , fc z ), and k = |k|. For simplicity, we have put 
72 = 73 in the original Luttinger Hamiltonian. In this 
Hamiltonian, a helicity defined byA = k- S//cisa good 
quantum number and can be used as a label for eigen- 
states. The helicity A can take values ±|, ±|; A = ±| 
correspond to the heavy- hole (HH) band and A = ±i 
to the light-hole (LH) band. Thus the energy EkA of 
the eigenstatestate with wavenumber k and helicity A is 
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expressed as 
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which only depends on k — |k|. By a unitary transforma- 
tion £/ k = e ies y e l< t>Sz ; w here 9, (ft are spherical coordinates 
of k, (k ■ S) is diagonalized as 

U k (k ■ S)UI = S Z = diag (| 1, -i -0 , (4) 

where fc = k/fc, and the Hamiltonian is diagonalized as 
well; 

U k H Ul = dmg{E kH , E kL , E kLl E kH ) . (5) 

Therefore, an eigenstate |kA) of the Hamiltonian H with 
hclicity A is expressed as 



M) = ul\\). 
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where |A) denote the eigenstate of S z with the eigenvalue 
A. 

We consider randomly distributed nonmagnetic im- 
puries with short-ranged potential: 



y(r) = ^^<5(r-R i ). 



(7) 



We shall calculate the Green function G = (z — -ff) -1 , 
H = Hq + V in the self-consistent Born approximation 
from the unperturbed Green function G' ' = (z — Hq) -1 . 
G (o) 

is diagonal in the helicity basis, with elements 



(z - E, 
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By taking the average (■ • • } 
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over the impurity distribution, the lowest-order term of 
the self-energy in the Born approximation is written as 
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where n = N/L 3 is a number density of the impurities. 

We shall show that this self-energy JSJ) is written as 
£<5kk"<5AA" > where E is a constant. By substituting © 
into ©, we get 

^(kA|VG (0 V|k"A")) Av 



<5 kk "^(A|C/ k ^|A'). 9 i?l,{A'|C/ k ,f/ k |A"}.(9) 
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We note that the Green function g k °,y is written as a k °) + 
b k °?(S 2 )yy, where ajj°> + \bf = (z - E k , L )~\ a™ + 



> k 

I tfe = ( z ~ -^fe'ff) -1 - Here ajj.7 and b k , depend only on 
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k' = |k'| and not on the direction of k'. The self-energy 
is then rewritten as 



(kA|VG (0 V|k"A") 



AV 



nV 2 
~W 

nV 2 



A" 



<5 kk » Y,\ X \ Uk K + • S) 2 Ul A"H 



Among various terms in (k' ■ S) 2 , the contributions from 
cross-terms such as k x k' y S x S v vanish by the summation 
over k . Likewise, in a term (k' x ) 2 S 2 , for example, (k' x ) 2 
can be replaced by \{k') 2 = g. Hence, (k ■ S) 2 can be 



replaced by \S 2 = §, and we get 
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By a similar argument, the full self-energy in the self- 
consistent Born approximation is seen to be diagonal in 
k and A: 



«kA|G|k"A")) AV 
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= 3fcA<5kk"<>AA", (12) 

5kff = S k ,§ =5k,-§> 3kL = Ski = 5k, ( 13 ) 
nV 2 



E = 
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where gfcA is a full Green function in the self-consistent 
Born approximation. The self-energy has several inter- 
esting properties; (i) it is diagonal in k and A, and (ii) it is 
independent of k and of A. It guarantees an assumption 
in the paper by Schliemann and Lossii. 

Let us calculate a vertex correction for a charge cur- 
rent vertex within the ladder approximation. The charge 
current operator J x is written as 



where { , } is the anticommutator. The lowest order term 
of the vertex correction for J x in the ladder approxima- 
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tion is given by 
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From eq. I|13(l . c/kA can be written as c/kA = a k + b k (S 2 )\\. 
This relation is useful in the calculation of the vertex 
correction as demonstrated below. 
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x J x {kx){a kl + b kl {k • S) 2 }[/t|A'). (17) 

Because the summand is an odd function of ki , this low- 
est order term of the vertex correction in the ladder ap- 
proximation vanishes identically. Therefore, higher or- 
der terms in the ladder approximation vanish as well, 
and thus we conclude that the vertex correction for the 
current J x in the ladder approximation vanishes. This 
cancellation is between the intermediate states with ki 
and — ki, i.e., due to parity. It is similar to the famil- 
iar example of vanishing vertex correction for the Fermi 
gas with spinless isotropic impurities. When we turn to 
the calculation of the spin Hall conductivity, we have to 
take into account the vertex correction in the charge cur- 
rent only. Therefore, because of this vanishing vertex 
correction, the spin Hall conductivity in the clean case 
reproduces the intrinsic value, which is obtained from a 
clean system without impurities from the outseliS* 1 ^* 3 ^. 
In other words, the result by taking a limit 1/r — > be- 
fore u) — > is identical with that by the reverse order of 
limits. The Luttinger model is free from a problem of 
order of limits. 

By inspection, we can generalize the above discus- 
sion; for any inversion-symmetric models with H(k) = 
H{— k), the vertex correction vanish identically for short- 
ranged scatterers. For example, even if we introduce the 
anisotropy in Luttinger parameters 72 7^ 73, the vertex 
correction remains zero. If the impurity potential is long- 
ranged, forward scattering is preferred and the vertex 
correction no longer vanishes. A simple calculation simi- 
lar to that for short-ranged scatterers shows that the re- 
sulting spin Hall effect is enhanced by a factor of r t s r H /r, 
where r tr H represents an effective transport lifetime for 
the spin Hall effect. 



This is to be contrasted with the Rashba Hamiltonian 
studied by Inoue et ah 24 ' 31 They studied the vertex cor- 
rection within the ladder approximation for randomly 
distributed nonmagnetic impurities with isotropic poten- 
tial. The result is remarkable; with the vertex correction, 
the spin Hall effect becomes zero in the clean limiljS 4 * 3 ^^, 
instead of e/8ir, the universal value without impurities 7 . 
In their paper, they calculated the charge-current ver- 
tex appearing in the correlation function between the 
charge current and the spin current in the Kubo formal- 
ism. In the clean limit, the vertex correction for the 
charge current in the ladder approximation is —1 times 
the spin-dependent part of the charge current operator; 
therefore, the total charge-current vertex becomes spin- 
independent, yielding a vanishing spin Hall conductivity. 
This vanishing result by Inoue et al. corresponds to the 
order of limits, to — > before 1/r — > 0, while the cal- 
culation of the intrinsic universal value by Sinova et al. 
corresponds to 1/t — > before to — > 0. The former and 
the latter limits correspond to L > / and L < I, respec- 
tively, where L is the system size and I is the mean free 
path. Nonetheless, it is not the end of the story. Two 
numerical calculations including disorder, one based on 
Kubo formula 34 and the other on scattering theory 2 ^, in- 
dicate that the spin Hall effect in the clean limit remains 
the universal value e/8n. It contradicts the analytical re- 
sult, and we do not have any answer to it. Here, we also 
note that Burkov et al22k formulated a theory of spin- 
charge coupled transport, and discuss that the spin Hall 
conductivity vanishes in a diffusive regime. Thus, a com- 
prehensive understanding of disorder effects for various 
strength of disorder is still to be desired. 

Though the Rashba model is different from the Lut- 
tinger model, our calculation on the Luttinger model has 
some implications on the debate on the disorder effect in 
the Rashba model. To interpret the vanishing spin Hall 
conductivity in the clean limit 2 ^, Inoue et al. discuss that 
diffuse scattering efficiently scrambles the precession of 
spins such that no net spin Hall current remains. This 
picture is too simplified, since the spin Hall conductivity 
is not necessarily zero for general systems with spin-orbit 
coupling. Indeed, in the Luttinger model the vertex cor- 
rection vanishes as we have seen. Furthermore, even for 
more general models with Rashba coupling, the spin Hall 
conductivity including the vertex correction is not nec- 
essarily zero in the clean limit. For example, instead of 
the simplest model of Rashba coupling 



h 2 k 2 

H = — h hX{k y o x - k x a v ), 
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let us take a model 

h 2 k 2 
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+ h(X + X 1 k 2 )(k y a x -k x cr y ) 1 (19) 



where A and Ai are constants. An extra term Ai is added 
here; this term should exist in general because it is al- 
lowed by symmetry. The vertex correction is calculated 
in the similar procedure as in refs. and we can 
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calculate the coefficient A' for the vertex correction. Af- 
ter lengthy but straightforward calculation, we can see 
that the spin Hall conductivity is nonzero in the clean 
limit, even if the vertex correction is included within 
the ladder approximation. Thus, for the simplest model 
of Rashba coupling, (|18J) . the complete cancellation of 
the spin Hall eSedM^ seems merely accidental, and 
not a consequence of any symmetries. This is also sup- 
ported by the result with long-ranged scatterers, where 
the spin Hall conductivity including the vertex correction 
no longer vanishes^. To summarize, we can say that the 
spin Hall effect is not necessarily suppressed to zero by 
the vertex correction in general. 



In conclusion, we consider an effect of spinless impuri- 
ties on the spin Hall effect in the Luttinger model. We 
calculated the vertex correction for the charge current 
within the ladder approximation. For short-ranged scat- 
ters, the vertex correction is zero, and the spin Hall effect 
reproduces the intrinsic value, calculated previously from 
the system without impurities. 
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